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PREFACE 


Thougli  this  rcporh  as  per  its  title,  deals  primarily  with  the  theoretieal  and  experimental 
investigation  of  acoustic  dyadic  sensors,  it  sets  the  groundwork  for  a  generalized  theory  of 
directional  acoustic  sensors.  Starting  with  the  Taylor  scries  of  the  acoustic  pressure,  it  moves 
quickly  from  the  basic  pressure  sensor  to  the  vector  sensor  and  then  on  to  the  dyadic  sensor. 

More  generalized  acoustic  sensors  arc  obtained  by  including  additional  terms  of  the  Taylor  scries. 
One  of  the  payoffs  of  more  advanced  sensors  is  improvement  in  the  bcamwidth.  With  a  vector 
sensor,  we  can  achieve  a  bcamwidth  of  105“  from  a  single  point.  Using  an  acoustic  dyadic 
sensor,  this  bcamwidth  can  be  decreased  to  65".  A  further  reduction  in  bcamwidth  is  anticipated 
as  the  order  of  the  acoustic  sensor  increases  (achieved  by  the  addition  of  more  Taylor  scries 
terms). 

This  report  shows  how  finitc-diffcrcncc  approximations  can  be  utilized  to  estimate  the  10  terms 
that  define  the  dyadic  sensor.  The  associated  realization  of  the  dyadic  sensor  takes  on  the  form  of 
19  judiciously  placed  pressure  sensors.  An  altcmativc  realization  using  18  accelerometers  and  I 
pressure  sensor  is  also  considered.  Both  realizations  involve  a  symmetrical  placement  of  sensors, 
which  allows  us  to  consider  the  sensors  as  collocated  with  a  common  phase  center  at  the  origin. 

A  detailed  analysis  is  carried  out  in  this  report  of  the  ciTor  that  arises  when  using  finite 
differences  in  an  environment  of  mismatched  sensors  that  have  nonzero  spacing  between  them. 
The  analysis  shows  that  it  is  not  prudent  to  employ  the  method  of  finite  differences  to  the 
estimation  of  second-order  partial  derivatives  of  the  pressure.  The  opposite  is  true  for  the 
estimation  of  first-order  partial  derivatives  (the  gradient  of  the  pressure).  For  example,  for  a 
+  0.2  dB  amplitude  mismatch  and  a  one-tenth  of  wavelength  spacing,  a  typical  mean  error  in 
estimating  a  first-order  partial  is  about  -0. 14  dB  with  a  standard  deviation  of  about  0.06  dB.  The 
error  analysis  leads  to  the  recommendation  that  an  acoustic  dyadic  sensor  be  realized  with  1 8 
accelerometers  and  a  pressure  sensor  at  the  origin.  This  will  eliminate  the  gradient  estimation 
error  and  significantly  reduce  the  Hessian  estimation  error,  since  the  Hessian  of  the  pressure  can 
now  be  estimated  by  first-order  finitc-diffcrcncc  approximations. 

This  report  also  shows  that  these  first-order  finitc-diffcrcncc  approximations  can  be  successfully 
performed  in  the  presence  of  quantization  and  wide  band  noise  effects.  A  worst-case  analysis 
shows  that  using  two  accelerometers  spaced  3  inches  apart,  one  can  achieve  a  normalized  RMS 
error  of  1 0%  when  attempting  a  first-order  finite  difference  of  a  50  Hz  plane  wave  in  a  2000  Hz 
band.  This  requires  a  signal-to-noisc  ratio  (SNR)  of  at  least  30  dB  and  assumes  that  the  noise  at 
each  accelerometer  is  uncorrclated.  This  can  be  done  using  a  1 6-bit  analog-to-digital  converter 
(ADC).  Further,  for  frequencies  up  to  2000  Hz,  the  required  SNR  drops  off  at  about  6  dB  per 
octave.  This  implies  that  the  same  3-inch  spaced  accelerometers  only  need  a  SNR  of  6  dB  to 
achieve  the  same  10  %  error  for  a  first-order  finite  difference  of  a  800Hz  plane  wave. 

A  partial  dyadic  sensor  (three  orthogonal  accelerometer  dipoles  and  a  pressure  sensor  at  the 
origin)  was  tested  at  Seneca  Lake  during  May  2001 .  The  experimental  in-water  beam  pattern 
measurements  confirmed  theoretical  predictions.  SITTEL  CORPORATION  also  demonstrated 
that  for  a  plane  wave  source,  dyadic  sensor  performance  can  be  achieved  by  advanced  signal 
processing  of  only  the  pressure  and  acceleration  measurements. 

This  report  was  prepared  for  James  McEachern  and  his  ONR  32 1 SS  team  under  ONR  contract 
number  N00014-01  -M-007.  The  Program  Officer  was  Jan  Lindberg.  This  is  the  Final  Report  and 
is  submitted  in  accordance  with  CDRL  AOOl,  Item  number  0001 AA-  Final  Report. 
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A  Theoretical  and  Experimental  Investigation 
Of  Acoustic  Dyadic  Sensors 


1.0  Introduction 

Current  research  in  the  area  of  directional  acoustic  sensors  has  shown  that  a  tri- 
axial  accelerometer  and  an  omnidirectional  hydrophone,  packaged  in  the  same  housing 
with  a  common  acoustic  phase  center,  can  provide  an  effective  directional  acoustic  sensor 
[  1 ,2,3].  Since  this  type  of  directional  acoustic  sensor  measures  both  the  scalar  acoustic 
pressure  and  the  x,  y  and  z  components  of  the  vector  part  of  the  acoustic  wave  field  (e.g., 
the  acoustic  particle  displacement,  velocity  or  acceleration),  it  has  been  called  an  acoustic 
vector  sensor  [4,5,6],  A  single  acoustic  vector  sensor  can  provide  up  to  6  dB  of  amay 
gain  against  isotropic  noise  relative  to  an  acoustic  scalar  sensor  or  omnidirectional 
hydrophone  [5].  This  vector  sensor  can  also  unambiguously  estimate  the  direction  of 
arrival  (DOA)  of  an  acoustic  source  [4].  Further,  this  vector  sensor  can  produce  a 
frequency-independent  beam  pattern  with  a  3-dB  beamwidth  of  105  degrees  [6,7], 

If  three  aceelerometers  and  a  pressure  sensor  (i.e.,  a  vector  sensor)  can  achieve  the  above 
results,  a  question  that  naturally  follows  is  how  many  more  individual  sensors  (e.g., 
accelerometers  and/or  pressure  sensors)  must  fit  into  a  given  sensor  housing  or  form 
factor  to  improve  upon  the  aforementioned  directionality  of  the  acoustic  vector  sensor? 

In  Section  2.0,  formulating  a  rigorous  mathematical  framework  for  the  general  theory  of 
directional  acoustic  sensors  provides  a  partial  answer.  From  this  theory  we  show  that  the 
next  logical  step  to  further  improve  vector  sensor  directionality  is  to  define  the  acoustic 
dyadic  sensor.  Section  3  shows  that,  for  a  general  linear  acoustic  wave,  19  judiciously 
placed  pressure  sensors  or  1 8  accelerometers  and  one  pressure  sensor  can  approximate 
the  acoustic  dyadic  sensor.  Section  4  consists  of  an  acoustic  analysis  of  a  theoretical 
dyadic  sensor,  and  Section  5  shows  how  an  acoustic  dyadic  sensor  can  be  used  to  create  a 
Multichannel  filter.  Experimental  results  from  a  Seneca  Lake  test,  conducted  in  May 
2001 ,  are  presented  in  Section  6.0.  Measured  beam  patterns  for  a  partial  dyadic  sensor, 
implemented  by  six  accelerometers  and  a  pressure  sensor,  are  presented  in  Section  6.1 . 

In  Section  6.2,  we  show  how  multi-channel  signal  processing  of  a  vector  sensor  can  be 
used  to  realize  the  performance  of  a  dyadie  sensor  when  the  pressure  field  is  assumed  to 
be  an  arbitrary  plane  wave.  We  also  provide  experimental  results  to  validate  this  claim. 
We  close  this  report  with  Section  7.0,  which  provides  a  summary  and  recommendations 
for  future  research. 

2.0  General  Theory  of  Directional  Acoustic  Sensors 

Let  us  consider  a  single  acoustic  sensor  located  at  some  measurement  point 
=  (,v„,y,),2„).  A  Taylor  series  for  the  scalar  acoustic  pressure  field  p{t,r)  about  this 

point  would  include  the  scalar  (tensor  of  order  zero)  pressure  pitJi)  as  the  zero-order 
term,  the  pressure  gradient/vector  (tensor  of  order  one)  Vp{t,r^^)  at  the  point  as  part  of  the 
first-order  term,  the  dyadic  (tensor  of  order  two)  V'Vp{t,rff)  at  the  point  as  part  of  the 


1 


second  order  temi,  and  so  on.  Using  this  Taylor  series,  we  define  a  general  class  of 
directional  acoustic  sensors  as  follows  [6].  A  scalar  acoustic  pressure  sensor  (e.g.,  a 
hydrophone)  will  be  referred  to  as  a  directional  acoustic  sensor  of  order  zero.  This  sensor 
only  measures  the  scalar  acoustic  pressure  at  the  point  /}, ;  its  Taylor  series  about 

assumes  that  the  acoustic  pressure  field  p{tj  )  about  that  point  is  independent  of  the 
field  point  f  =  (x,  v,^).  An  acoustic  vector  sensor  will  be  refeiTed  to  as  a  directional 
acoustic  sensor  of  order  one.  This  sensor  measures  both  p{tj\^)  and  the  pressure 
gradient/vector  Vp(t,r„)  at  the  point  r„ ;  its  Taylor  series  about  r„  assumes  that  the 
acoustic  pressure  field  p(t,r  )  about  that  point  is  a  linear  function  of  the  field  point  f  . 
Similarly,  an  acoustic  dyadic  sensor  will  be  referred  to  as  a  directional  acoustic  sensor 
of  order  two.  This  sensor  measures  p(t,r,,),  Vp(t,/}J  and  the  dyadic  VV/;(/,r^)  at  the 

point  r„ ;  its  Taylor  series  about  /j,  assumes  that  the  acoustic  pressure  field  p(t,r)  about 
that  point  is  a  quadratic  function  of  the  field  point  r  . 

D’ Spain  [8]  was  the  first  to  point  out  a  qualitative  relationship  between  the  Taylor 
series  expansion  of  the  pressure  field  and  a  vector  sensor.  Following  his  observation,  we 
consider  a  single  acoustic  sensor  located  at  some  measurement  point  r„  =  (x„,  jn,Zo).  The 
Taylor  series  for  the  pressure  p(x,y,z,t)  expanded  about  the  point  (x,, ,  )  is 


V,  z,  t)  — /?(x()  ,yf),.2:o,r)+  (a  “*"0’  “o)-,  ’  To  >  ’  0 

ox  o  V  o  z 

+  -  -  ^'n  }-^  +  (y-  Vo  +  (^  -  ^0  p(^'0 .  To  Ao ,  0  +  ^3 


where  y>(xn,y(,,z„,r)  following  the  brackets  means  that  the  partials  within  the  brackets 
are  to  operate  upon  p  at  the  point  (x(,,yo’-o )’  the  brackets  indicate  that  the 
expansion  of  the  quantity  within  is  to  be  a  trinomial  expansion  except  that 
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is  to  be  replaced  by 


and  terms  like 
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are  to  be  replaced  by  1 .  Let  us  define  the  column  vectors 
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and  the  Hessian  matrix 
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Using  (5)  through  (7)  in  ( 1 )  and  dropping  the  remainder  term  7?, ,  we  obtain 

p{x, y, zj)=  p{x, ,  Vo ,z,„t)  +  {r-)p  Vp{x„ , j) 

+  T  t  -  fii )  ’  .17.  ^  ^0  ^  -  !ii) 

The  function  ^(x,  v,5,t)  is  the  pressure  estimate  generated  by  the  dyadic  sensor  and  is 
capable  of  extrapolating  the  acoustic  pressure  field  beyond  the  measurement  point  so 
that  it  actually  knows  this  field  at  every  point  inside  a  sphere  of  radius  7?  =  |r  -  rj,! ,  where 

R  is  defined  by  the  type  of  directional  sensor  and  the  error  associated  with  the 
extrapolation.  The  Taylor  series  (8)  is  the  fonnula  used  to  do  the  wave  field  extrapolation 
with  some  specified  error  s(f,r).  For  example,  a  directional  sensor  of  order  zero,  can 
only  measure  the  pressure  at  F„ ,  so  its  estimate  of  the  field  beyond  this  point  is 
p(/,r)=  p(/,rf,)  and  the  corresponding  estimation  or  extrapolation  error  is 
£•(/,?)  =  p(/,F)-  p(t,r,i).  If  the  error  is  required  to  be  small  (less  than  10%),  then  R  will 

be  small.  This  implies  that  the  aperture  (2R)  of  a  single  pressure  sensor  is  small,  so  by  the 
theory  of  spatial  Fourier  transfoirns  [9],  the  zero-order  sensor  is  essentially 
omnidirectional.  However,  a  sensor  of  order  one  (vector  sensor)  measures  both  the 
pressure  and  pressure  gradient  at  /j, ,  so  its  estimate  of  the  field  beyond  this  point  is 

p{uF)  -  p{i,F\^)  +  (F  -  F\^)-Vp{t,F\^)  and  the  corresponding  estimation  error  is 
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£{t,r)=  p{t,r)-  p{^,r  ) 
e{t,  r )  =  p{t,  r )  -  [p{t,  Fo )  +  (F  -  Fq  )  •  Vp{t,  )] 


(9) 


For  the  same  error,  the  vector  sensor  will  have  a  larger  aperture  than  the  scalar  sensor,  so 
it  should  be  more  directional.  Thus,  the  dyadic  sensor  should  be  better  than  the  scalar  and 
vector  sensors  at  estimating  the  field,  so  it  should  be  even  more  directional. 


toward  the  origin  of  a  rectangular  coordinate  system. 

For  the  case  of  a  plane  wave  pressure  field  (refer  to  Figure  1),  we  define  the  mean- 
squared  error  in  estimating  the  pressure  field  beyond  F^,  by 


j  IT  T 

MSE  =  —\\\s{t,7fdtdp,  (10) 

where  T  is  a  suitable  integration  time  (e.g.,  the  temporal  period  of  a  plane  wave)  and  p  is 
the  angle  between  h  and  F  (a  unit  vector  in  the  direction  of  F ). 
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Figure  2.  Normalized  mean-squared  estimation  error  (equation  (19))  vs.  R! X  . 
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Dividing  (10)  by  the  average  power  in  the  acoustic  pressure  field  we  obtain  the 
normalized  MSE.  Figure  2  shows  the  normalized  mean-squared  error  as  a  function  of 
R! Pi  for  the  scalar,  vector  and  dyadic  sensors.  For  a  specified  normalized  MSE  of  10%, 
notice  that  the  acoustic  aperture  (2R)  of  the  scalar  sensor  is  about  A /1 0,  whereas  the 
apertures  for  the  vector  and  dyadic  sensors  are  X/Zmd  XU  ,  respectively.  Figure  3  is 
scaled  to  graphically  depict  the  situation. 


I  Taylor  Series  Extrapolation  of  Pressure  Field  =  Volumetric  Spherical  Sensors 


scalar 


dyadic 


Figure  3.  Spherical  coverage  for  10%  error. 


We  are  now  in  the  position  to  answer  the  question  posed  in  the  Introduction:  If  three 
accelerometers  and  a  pressure  sensor  (i.e.,  a  vector  sensor)  can  achieve  the  above 
results,  a  question  that  naturally  follows  is  how  many  more  individual  sensors  (e.g., 
accelerometers  and/or  pressure  sensors)  must  fit  into  a  given  sensor  housing  or  form 
factor  to  improve  upon  the  aforementioned  directionality  of  the  aeoustic  veetor  sensor? 
Recall  that  the  dyadic  sensor  measures  p{t,r^),  V/i(?,ro)  and  the  dyadic  VV;?(t,r„)  at  the 
point  To .  According  to  (6),  (7),  and  (8),  we  must  measure  13  terms  involving  the  pressure 

and  the  first  and  seeond  partial  derivatives  of  the  pressure.  Due  to  the  symmetry  of  the 
mixed  partials  in  (7),  the  number  of  distinct  terms  reduces  to  10.  The  next  section 
demonstrates  that  these  10  ternis  can  be  estimated  by  19  judiciously  spaced  pressure 
sensors,  or  1 8  judiciously  spaced  aceelerometers  and  1  pressure  sensor  placed  at  the 
origin. 

3.0  Realization  of  an  Acoustic  Dyadic  Sensor 

The  goal  of  this  section  is  to  show  how  finite  difference  approximations  can  be  used 
to  estimate  the  1 0  tenus  that  define  the  dyadic  sensor;  that  is,  the  scalar  pressure,  the  3 
terms  of  the  pressure  gradient,  and  the  6  distinct  terms  in  the  Hessian  matrix.  We  will 
first  show  how  to  realize  the  dyadic  sensor  by  using  19  pressure  sensors,  then  we  will 
provide  a  dyadic  sensor  realization  that  uses  1 8  judiciously  spaced  accelerometers  and 
one  pressure  sensor  at  the  origin.  Both  realizations  utilize  a  symmetrical  placement  of 
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sensors,  which  allows  us  to  consider  the  sensors  as  collocated  with  a  common  phase 
center  at  the  origin. 

3.1  Finite-Difference  Approximation  of  the  Pressure  Gradient 

Figure  4  illustrates  an  arrangement  of  pressure  sensors  useful  for  estimating  the 

first  partial  derivatives  — ,  — and  —  at  (0,0,0).  The  finite-difference  approximations 

5  v  cz 

of  the  pressure  gradient  are  given  by 


ep(r,0.0,0)  ^  0,0) ,0,0)]/ 2/1 

dx 

(11a) 

o/7(r,0,0,0)  ^  o,/7,0)-/?(t,0,-/i,0)]/2/7 

dy 

(11b) 

dp(f, 0,0,0)  ^  o,0,/7)  p(t,0, 0, -/?)]/ 2/i 

dz 

(11c) 

The  errors  of  the  above  approximations  are  on  the  order  of  h' .  Note  that  6  pressure 
sensors  are  required  to  realize  the  pressure  gradient. 

Figure  4.  The  arrangement  of  six  pressure  sensors  for  the 
finite-difference  approximation  of  the  pressure  gradient. 
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3.2  Finite-Difference  Approximation  of  the  Hessian  of  the  Pressure 


Figure  5  illustrates  an  arrangement  of  pressure  sensors  useful  for  estimating  the 

second  partial  derivatives  ^ ^  and  at  (0,0,0).  The  finite-difference 

(2.V'  rv" 

approximations  of  these  partials  arc  given  by 


a-/;(0A0,t)  ^  pjliAOj)-  2 p(0,0,0, i)+p{-  /7,0,0, / ) 


<9V(0,0,0,/)  ^  /;{0,/?,0,/)-2;?(0,0,Ot)+  p{0-h,0,t) 

2  ^  j  2 

ov  h 

d' p{0fi.0j)  _  ;;(0,0,/?,t)-2p(0,0,0,?)+  p{0,0-h,t) 


dz" 


Jr 


(12a) 

(12b) 

(12c) 


The  errors  of  the  above  approximations  are  on  the  order  of  Jr .  Note  from  (12)  and  from 
Figure  5  that  there  are  7  pressure  sensors  involved  in  the  estimation  of 

^  and  — ^  at  (0,0,0).  The  pressure  sensor  at  the  origin  is  also  utilized  to 

dv  &■ 

measure  the  acoustic  pressure  at  (0,0,0).  This  pressure  is  the  first  term  in  the  pressure 
estimate  given  by  (8)  where  =  v,,  =  s,,  =  0  .  The  other  six  sensors  in  Figure  5  are  the 
same  sensors  that  provide  data  for  the  estimates  appearing  in  Equations  (11). 


Figure  5.  Sensor  placement  for  the  estimation  of  the  diagonal 
elements  of  the  Hessian  of  the  pressure  at  the  origin. 
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Figure  6  illustrates  an  arrangement  of  pressure  sensors  useful  for  estimating  the 
mixed  second  partial  derivatives 


^  ^  ^  and-^-^  at  (0,0,0).  The  four  “star’' 


dxdv  dxdz 


ovcz 


sensors  can  be  used  to  estimate  d' 


,  the  four  “diamond”  sensors  can  be  used 


to  estimate  d~  pjdxdz  ,  and  the  four  “square”  sensors  can  be  used  to  estimate 


)-plc 


cF  p  dvdz 


.v„v=r--=0 


.  The  finite-difference  approximations  of  these  partials  are  given  by 


dxdy 

d-p{0,0,0j)  .  pfeA-rj) 

\-p\ 

(4,0 -T’O 
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'+  p{ 
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(13a) 


(13b) 


(13c) 


The  errors  of  the  above  approximations  are  on  the  order  of  Jr .  Note  from  (13)  and  from 
Figure  6  that  there  are  a  total  of  12  pressure  sensors  involved  in  the  estimation  of 

2  2  ^2 

Z—E-  ^  ZJi.  and  — ^  at  (0,0,0).  This  makes  a  total  of  19  sensors:  7  for  the  estimation  of 
dxdy  dxdz  dydz 

the  gradient  and  the  diagonal  terms  of  the  Hessian  and  12  for  the  estimation  of  the  off- 
diagonal  terms  of  the  Hessian. 


d-p 

dxd}’ 
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Figure  6.  Sensor  placement  for  the  estimation  of  the  off-(iiagonal 
elements  of  the  Hessian  of  the  pressure  at  the  origin. 
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3.3  Realization  of  a  Dyadic  Sensor  by  Means  of  Accelerometers 


From  Newton’s  Second  Law  of  Motion,  the  pressure  gradient  and  acoustic  particle 
acceleration  vector  are  related  by 

Vp  =  -pa  (14) 


Taking  the  vector  gradient  of  both  sides  of  (14)  results  in  the  dyadic/Hessian 


where 


Wp  -  -V(yC>  r/)  =  -pVa  -aVp»  -pVa 


(15) 
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dx 

ox 
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dx 
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ocy. 
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az 

dz 

(16) 


is  the  gradient  of  ^ ,  /?  is  the  mass  density  and  Vp  is  the  gradient  of  the  mass  density  at 
the  origin.  Note  the  approximation  appearing  in  (15).  According  to  Camp  [10],  a 
reasonable  acoustic  approximation  is  to  neglect  Vp  at  the  origin.  Equations  (14)  and 

(15)  arc  two  of  the  three  quantities  measured  by  a  dyadic  sensor  (pressure  being  the  third 
one).  First-order  finite  difference  approximations  can  be  used  to  estimate  the  elements  of 

(16) .  This  would  require  18  accelerometers  symmetrically  arranged  as  in  Figures  5  and 
6.  Like  the  pressure  sensor  realization  of  the  dyadic  sensor,  the  accelerometer  realization 
utilizes  a  symmetrical  placement  of  accelerometers,  which  allows  us  to  consider  the 
accelerometers  as  collocated  with  a  common  phase  center  at  the  origin. 


3.4  Comparison  of  Pressure  and  Accelerometer  Realizations  of  the  Dyadic  Sensor 

Table  1  provides  comparisons  between  the  pressure  and  acceleration  methods  of  dyadic 
sensor  realization: 


Table  1 


Pressure  Realization 

Accelerometer  Realization 

Requires  1 9  pressure  sensors  with  one  at  the 
origin.  See  Figures  5  and  6  for  placement. 

Requires  19  sensors  (18  accelerometers  and 
one  pressure  sensor).  Sec  Figures  5  and  6  for 
placement. 

One  hydrophone  at  origin  measures  pressure 

One  hydrophone  at  origin  measures  pressure 

Estimates  the  pressure  gradient  by  first-order 
finite  differences. 

Measures  the  pressure  gradient.  Sec  Eq.  (14). 

Docs  not  require  knowledge  of  p  at  the 
measurement  point  to  estimate  the  Hessian. 

Requires  knowledge  of  p  at  the  measurement 
point  to  estimate  the  Hessian. 

Requires  first-order  and  second -order  finitc- 
diffcrcncc  approximations 

Requires  only  first-order  finite-difference 
approximations 
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4.0  Analysis  of  a  Theoretical  Acoustic  Dyadic  Sensor 

This  section  will  examine  the  effects  of  finite-difference  approximations  on  the 
estimation  of  the  pressure  gradient  and  Hessian  (realization  by  pressure  sensors).  It  will 
also  examine  the  effects  of  finite-difference  approximations  on  the  estimation  of  Va 
(realization  by  accelerometers). 

4.1  Pressure  Gradient  Estimation 


In  the  undewater  and  in-air  communities,  the  pressure  gradient/vector,  needed  for 
a  vector  sensor,  is  found  in  a  variety  of  different  ways.  As  discussed  above,  one  can 
measure  the  acoustic  particle  acceleration  vector  directly  and  use  (14)  and 


knowledge  of  the  mass  density  at  the  measurement  point  Ff,  to  obtain  a  fairly  accurate 
value  of  the  pressure  gradient  Vp{t,r^,) .  Alternatively,  one  could  measure  the  acoustic 
particle  velocity  vector  v(t,/;,)  directly,  perform  a  time  derivative  to  get  the  acceleration 

vector  and  proceed  as  above  to  get  the  pressure  gradient.  (The  accuracy  of  this  method 
would  depend  on  how  well  one  implemented  the  time  derivative).  If  one  decides  to  use 
scalar  pressure  sensors  to  approximate  the  pressure  gradient,  then  the  finite  difference 
approximations  in  ( 1 1 )  would  have  to  be  computed.  However,  great  care  must  be 
exercised  when  approximating  the  pressure  gradient  by  (1 1).  For  example,  (11a) 
approximates  the  x-component  of  the  pressure  gradient  at  the  measurement  point 

=  (0,0,0) .  (Refer  to  Figure  4).  If  the  pressure  field  was  the  arbitrary  (narrowband  or 


broadband)  plane  wave 


p{tJ-)  =  f 


^  fi-r^ 
t  + - 


V  c  J 


(17) 


which  is  assumed  to  be  propagating  toward  the  origin  of  the  coordinate  system  in  Figure 
1 ,  then  the  exact  value  of  the  temporal  Fourier  transform  of  the  pressure  gradient  would 
be 


VP[co,r)=  jF[co)ke'^  '  =  jP{co,r)k  .  (18) 

Here,  the  wavenumber  or  propagation  vector  k  =  kh  has  the  magnitude  k  =  col c  =  2n I X 
and  the  direction  n  ,  where  h  is  the  unit  vector 

h  =  sin  (pj  cos  BjX  +  sin  (f)j  sin  Oj-y  +  cos  ^.,  z .  (19) 

The  position  vector 

r  =  xx  +  yy  +  zz  =  /-sin^cos^ -i- rsin^cos^  +  rcos^zis  (20) 
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is  the  vector  from  the  origin  to  an  arbitrary  field  point  r  —  (x,  v,^:)  in  rectangular 
coordinates  or  r  =  (/•sin^cos(9,/'sin^sin6?,/-cos^)  in  spherical  coordinates.  The 
temporal  Fourier  transform  of  the  plane  wave  pressure  field  (17)  is  found  by  using 


(21) 


Now  if  we  use  (1  la)  to  approximate  the  x-component  of  the  pressure  gradient  at  the 
origin  (the  point  =  (0,0,0) ),  then  the  temporal  Fourier  transform  of  this  approximation 
can  be  expressed  as 


dP{co,r,.)  ^  /  _  X  -  /  _  .  /  ssinf/r  h) 


ax 


(22) 


Since  the  exact  value  of  is  given  by 


P^icoJ-,)=jk^F{co),  (23) 

as  /i  — >•  0,  the  estimate  given  by  (22)  approaches  the  exact  value  given  by  (23).  A  good 
measure  of  the  quality  of  the  gradient  estimation  is  the  decibel  error  ,  namely, 


f  =20  log 


(24) 


Since  h  =  ,  where  Ar  is  the  sensor  spacing,  we  have  from  (22)  and  (23)  that 

^  =201og|sinc  (A',,Aa-/2J  (25) 


Figure  7  illustrates  the  error  f  as  a  function  of  Aji'/T  for  the  case  =  90°  and  O-i  =  0° . 
An  eiTor  of -0.5  dB  will  occur  for  a  spacing  of  Ax  =  0.186/1  or  0.186  wavelengths.  A 
spacing  of  one-tenth  of  a  wavelength  will  yield  an  error  of -0. 14  dB.  Another  popular 
measure  of  the  quality  of  estimation  is  the  percent  error  defined  as 


S.  = 


xlOO 


(26) 


The  relationship  between  and  f  is  illustrated  in  Figure  8.  For  example,  a  value  for  # 
of  0.4  dB  corresponds  to  a  of  about  4.7%. 
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Figure  1.  The  error  /  in  estimating  the  partial  derivative  of  the  pressure  with 
respect  to  x  (at  the  origin).  (/)j  =  90  and  6j  —  0^ . 


Figure  8.  Relationship  between  percent  en*or  and  dB  error. 


Notice  that  (25)  is  frequency  dependent  because  of  the  presence  of  and  is  valid  only 
when  the  two  hydrophones  in  the  x-axis  dipole  (Refer  to  Figure  4)  are  exactly  matched  in 
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their  amplitude  and  phase  frequency  responses.  In  practice,  this  is  seldom  the  case,  so  we 
must  modify  our  analysis  to  include  the  effects  of  mismatch  in  the  frequency  responses  of 
the  hydrophone  and  corresponding  signal  conditioning  associated  with  each  hydrophone 
in  the  dipole.  The  temporal  Fourier  transform  of  ( 1 7)  is 

P{coJ-)  =  F{co)e‘''''^  (27) 


However,  when  the  pressure  waveform  impinges  on  the  sensor,  the  sensor  would 
transform  this  pressure  to  the  analog  voltage  Fourier  transform 

V[(o,?)  =  K{co)p{coJ-)  (28) 


The  quantity  K(rt))  is  the  composite  (sensor  sensitivity  plus  electronics)  frequency 
response  of  the  sensor.  Taking  the  temporal  Fourier  transform  of  (I  la)  and  incoiporating 
the  composite  frequency  responses  of  the  two  sensors,  we  arrive  at 


V{coJ-,)=^F{co) 


K,(ry)e''''' 


2/7 


(29) 


Let  the  mismatch  between  the  two  sensors  be  described  by  the  quantity 

M(ru)  =  ^^^  =  /7(ru)exp[/^(ry)],  (30) 

K,(w) 

where  //(ru)  and  i^{(o)  are  the  amplitude  and  phase  mismatch,  respectively. 


Utilizing  (30)  in  (29)  gives 

K|(7y)  ^  ’  2/7 


(31) 


The  magnitude  of  (3 1 )  is 

HoJ-A 


kM 


■yj\  -  2 /j cosilk  J't-(^)  +  /d' 


(32) 


The  dB  gradient  estimation  error  has  the  value 

•^1  -2/ycos(2A',  /7  -4')+  p' 


f  -  20  log 


Ikh 


(33) 
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When  //  =  1  and  =  0 ,  there  is  no  amplitude  and  phase  mismatch  and  (33)  reduces  to 
(25).  The  error  f  becomes  infinite  when  k\  -  0 .  However,  A',  =  0  means  that  the 

pressure  waveform  is  not  a  function  of  x ,  and  we  would  not  be  estimating  the  partial 
derivative  of  the  pressure  with  respect  to  a'  since  it  is  zero.  From  Equation  (1  la)  we 
observe  that  the  separation  between  the  two  sensors  along  the  x  axis  is  Aa'  =  2h  . 

Assuming  that  -  0  (no  phase  mismatch)  and  noting  that  A  ^  =  ImiJX,  we  can  write  (33) 
as 


6  =  20  log 


.yl-2//cos(2;r(/?,.AA'//l))  +  //” 
2;r(/7^,AA7/l) 


(34) 


Figure  9  is  a  contour  plot  of  f  versus  Ax//l  and  //  for  the  case  =  sin  (f)-,  cosOj-  =  ±1 
(f'  is  an  even  function  of  ).  The  contours  shown  con'espond  to  ft'  =  ±1,  ±  0.5  and  0 


Figure  9.  Contour  plot  of  the  error  ft  as  a  function  of  p  and  AxA  for  the 
case  =  +1 . 


dB.  From  the  figure  we  see  that  to  achieve  a  gradient  estimation  error  of -0.5  dB  (-5.6%) 
when  the  sensor  amplitude  mismatch  is//  =  -1  dB,  we  need  a  sensor  spacing  of  0.1 
wavelengths.  This  is  46%  smaller  than  the  required  relative  spacing  with  no  mismatch. 
Figure  10  is  a  contour  plot  for  the  case  n^.  =  ±4 .  Note  that  for  the  same  S  and  //  the 

sensor  spacing  is  now  0.2  wavelengths  or  double  the  previous  value.  For  an  arbitrary  , 
the  new  spacing  would  be  0.l/|7/,,|  wavelengths.  Note  that  the  spacing  increases  as  [n  J 
decreases.  Hence,  j/i  J  =  1  is  the  most  restrictive  case  and  Figure  9  should  be  used  to 
detemiine  the  required  sensor  spacing. 
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Figure  10.  A  contour  plot  of  the  crroiV  for  the  ease 


4.2  Estimation  of  the  Pure  Second-Order  Partial  Derivatives  of  Pressure 


This  section  performs  an  analysis  of  the  error  that  results  in  utilizing  finite-difference 
approximations  to  estimate  the  pure  second-order  derivatives  of  aeoustic  pressure.  The 
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•  - ► 


Ks  K,  K, 
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Figure  1 1 .  Arrangement  of  pressure  sensors  for  the  finite-difference 
approximation  of  the  pure  second-order  partial  pj  OX^  at  the  origin. 


errors  that  will  be  considered  are  those  due  to  non-zero  sensor  spacing  and  to  mismatch 
in  the  frequency  responses  of  the  various  sensors.  We  will  be  concerned  with  the  second- 
order  partial  derivative  d'  p{0,0,0,t)/dx'  .  The  results  of  error  analyses  of  the  two  other 
pure  second-order  partial  derivatives  would  be  analogous.  Figure  1 1  shows  the 
arrangement  of  the  pressure  sensors. Starting  with  the  estimation  equation  (13a)  we  would 
arrive  at  an  equation  analogous  to  (29),  namely, 
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(35) 


K,  exp(,/A-^/?)-  2K.  +  K,  exp(-  jkji) 


If  K-^  =  K2  =  (no  mismatch),  then  (35)  reduces  to 

V^v(®)  =  “/',’'^i^(^'^)siric‘(A'v  /?/2)  (36) 

From  Figure  5  we  observe  that  the  relative  sensor  spacing  is  Aa'  =  h  .  The  error  <5  is 
given  by 


^  =  20  log 


d'  pi^,o))ldxdy 
d'  P^,co)l dxdy 


40  log|sinc(A, AY/2)j 


(37) 


Equation  (37)  approaches  zero  as  Av  approaches  zero.  A  comparison  of  (37)  with  (25) 
reveals  that  for  a  given  sensor  separation  Ay  ,  the  en'or  in  estimating  the  pure  second- 
order  partial  derivative  is  twice  the  error  in  estimating  the  first-order  partial  derivative. 
The  problem  is  more  complicated  when  the  sensors  are  mismatched.  The  magnitude  of 
(35)  has  the  fonn 

=  (38) 


where 


r  = 


1 

■  1  ' 

A, /A, 

P2 

1  -2cos(A,.A)  cos{2kJi) 

-  2  cos(A  ^  A )  4  -  2  cos(A  ^  h ) 

cos{2kyh)  -2cos(A^/?)  1 


For  this  more  general  case,  (37)  becomes 


(f  =  20  log 


(39) 


(40) 


(41) 


If  each  element  in  /j  is  one,  then  (41)  reduces  to  (37).  In  the  current  error  analysis  we 

have  ignored  the  mismatch  in  phase  response  between  sensors.  A  more  complete 
analysis  would  result  in  addition  of  a  phase  terms  in  the  arguments  of  the  cosines  that 
appear  in  (40).  To  see  the  effect  of  sensor  mismatch  on  the  dB  error  ,  let  us  begin  by 
computing  the  relative  sensor  spacing  to  achieve  an  S  of -0.4  dB  (4.7%  error)  when 
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(/>i-  =  90  and  O,  =  0' .  ]f  we  use  (37)  we  will  obtain  A,v//l  =  0.1 18  .  Figure  12  shows  a 
contour  plot  of  ft  in  dB  versus  the  mismatch  parameters  //^  and  //,  for  the  special  case 
(/),  =  90  , 0,  =  0  and  Av/ 2  =  0. 1 1 8 .  Observe  from  Figure  1 2  that  when  //,  =  //,  =  0  dB , 

the  value  of  /  is  -0.4  dB  as  expected.  We  see  from  the  figure,  that  sensor  mismatch  may 
improve  the  error  (note  the  0  dB  error  contour).  However,  the  error  can  get  much  larger. 
For  example,  if  //,  =  -0.5  and  //,  =  +0.5  dB,  the  value  of  f  increases  to  about  2  dB. 


Figure  12.  Contour  plot  for  as  a  function  of  the  mismatch  parameters 
M:  and  p;,  for  the  special  case  of  (|)t  =  90'\  Gy  =  C'  and  Ax/^=  0.118. 


4.3  Estimation  of  the  Mixed  Second-Order  Partial  Derivatives  of  Pressure 


This  section  performs  an  analysis  of  the  error  that  results  in  utilizing  finite-difference 
approximations  to  estimate  the  mixed  second-order  derivatives  of  acoustie  pressure.  The 
errors  that  will  be  considered  are  those  due  to  non-zero  sensor  spacing  and  to  mismatch 
in  the  frequency  responses  of  the  various  sensors.  We  will  be  concerned  with  the  second- 
order  partial  derivative  d'  p{0fi,Q,t)ldxdy  .  The  results  of  error  analyses  of  the  other 
mixed  second-order  partial  derivatives  would  be  analogous.  Figure  13  shows  the 
arrangement  of  the  pressure  sensors. 
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Figure  13.  Arrangement  of  pressure  sensors  for  the  finite-difference 
approximation  of  the  mixed  second-order  partial  d"  pjdxdy  at  the  origin. 


Starting  with  the  estimation  equation  (13a)  we  would  arrive  at  an  equation  analogous  to 
(35),  namely, 


.  V  ^  A-K~k.W 

^7  i  \  T^i  -A,e  +K^e  -Jx^e 

V,,  \co)  =  F[co) - - —  (42) 

h 


If  =  A,  =  A,  =  A,  (no  mismatch),  then  (42)  reduces  to 

Vv,.(ry)  =  A|F(<y)sinc(^^.  /2/2)sinc(A',.  /i/2)  (43) 


The  error  A  takes  on  the  forrn 

A  =  201og|sinc(A',,/7/2)sinc(A',,/2/2]| 


(44) 


approaches  zero  as  h  approaches  zero.  The  problem  is  more  complicated  when  the 
sensors  are  mismatched.  The  magnitude  of  (42)  has  the  form 

where 
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(46) 


1 

■  1  ■ 

K,JK, 

A: 

K,/K, 

/■h 

_KJK,_ 

•cos 


I  -COs(/i,/7) 

■(*,*)  I 

-  cos{k  Ji )  cos((A'  ,  -  k )h  ) 

cos  ((*,  +*,>')  -C0s(A',/7) 


-  cos{kJi) 
cos((a,  -  a,  )?) 

] 

-  COs(A,./7) 


cos((a^+A,,>} 

-COs(A^/7) 

-cos(a,  77) 


(47) 


The  error  f  becomes 


=  201og| 


(48) 


When  the  vector  //  contains  all  ones,  (48)  simplifies  to  (44).  In  oi'der  to  gain  some 
understanding  of  the  behavior  of  (48),  let  us  define 


f 

dB(/)s201og 

V 


K 


I  J 


i  =  2,3,4 


(49) 


We  will  choose  a  positive  number  dBmax  and  declare  eaeh  dB(/)  to  be  a  random  variable 
uniformly  distributed  on  the  interval  (-dB„,ax,+dB„,ax)-  Onee  sample  values  are  obtained, 
the  sensor  sensitivities  are  given  by 


j'dli(/)\ 

K,  =K,W  >,  7  =  2,3,4  (50) 

Figure  14  shows  a  plot  of  the  error  (48)  versus  the  sensor  separation  in  terms  of  a  fraction 
of  a  wavelength  for  dB,„ax  going  from  0.1  to  2  dB  in  steps  of  0. 1  dB.  The  values  of 
(f)j  and  9j  are  90"  and  45",  respectively.  Equations  (49)  and  (50)  were  applied  1 00  times 
for  each  dBmax  and  the  maximum  error  plotted.  Figure  14  reveals  some  important 
features  of  the  estimation  error.  As  the  relative  sensor  separation  hj X  becomes  larger, 
the  error  is  less  sensitive  to  the  degree  of  sensor  mismatch.  Note  that  at  half-wavelength 
spacing  the  band  of  en  or  between  dB,„ax  =  0.1  dB  and  dB^ax  =  2  dB  is  narrow.  Flowever 
for  hj X  <  0.2 ,  the  relative  error  is  highly  sensitive  to  sensor  mismateh,  the  error  reaching 
unsatisfactory  values  very  quickly.  In  summary,  for  small  A/T  the  error  is  sensor- 
mismatch  limited  and  for  large  hj X  the  error  is  sensor-spacing  limited.  It  appears  from 
Figure  14  that  the  sensor  spacing  should  be  somewhere  in  midrange.  Figure  15  gives  a 
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4.4  Summary 

This  section  has  analyzed  some  of  the  errors  that  occur  when  finite-differences  are  used 
to  estimate  the  first-order  and  second-order  partial  derivatives  of  the  acoustic  pressure. 
The  errors  considered  are  those  caused  by  a  ( 1 )  finite,  nonzero  spacing  Ar  between 
sensors  and  (2)  an  amplitude  mismatch  among  the  composite  frequency  responses  of  the 
various  sensors.  Expressions  for  these  errors  were  derived  and  plots  were  generated  that 
illustrated  the  dependency  of  error-behavior  on  relative  sensor  spacing  and  dB-mismatch. 

This  section  ends  by  describing  pictorially  the  statistical  behavior  of  the  estimation 
error.  Let  (j)j  =  90  and  6j  =  45  .  Let  us  first  focus  on  estimating  d' p[0fifij)ldxdy . 
Figure  16  shows  a  histogram  of  the  estimation  error  (equation  (48))  at  a  sensor  spacing  of 
0. 1  wavelengths.  All  of  the  sensors  were  assumed  to  be  within  ±  0.2  dB  of  each  other  (a 
uniform  probability  distribution  for  each  was  assumed).  The  histogram  represents  a 

compilation  of  2000  sample  realizations.  Also  illustrated  in  Figure  16  is  the  mean  value 
and  standard  deviation  of  the  estimation  error  for  sensor  spacing  from  0. 1  to  0.5 
wavelengths.  At  a  spacing  of  0. 1  wavelengths,  the  mean  emor  is  about  0.25  dB. 

Flowever,  the  standard  deviation  is  1  dB.  Note  that  as  the  spacing  increases,  the  mean 
error  increases,  whereas  the  standard  deviation  decreases. 

Figure  17  contrasts  the  estimation  of  the  first-order  partial  derivative  (5j!7(0,0,0,/)/(5x  and 
the  second-order  partial  derivative V(0,0,0,r)/dv'  .  For  both  cases,  =  90"  and 
6j  =  0  .  Observe  the  plots  in  the  upper  portion  of  Figure  1 7.  The  mean  error  (in 
estimating  the  second-order  partial,  refer  to  equation  (41))  at  the  right  is  twice  the  mean 
error  (in  estimating  the  first-order  partial)  at  the  left.  The  lower  portion  of  Figure  17 
shows  the  standard  deviation  of  both  cases.  In  the  first-order  partial  case,  the  standard 
deviation  is  constant  with  respect  to  sensor  spacing  and  has  a  small  value  of  0.0571 .  In 
the  second-order  partial  case,  the  standard  deviation  decreases  exponentially  from  0.66  at 
a  spacing  of  0. 1  wavelengths  to  0.065  at  a  spacing  of  0.5  wavelengths.  At  a  sensor 
spacing  of  0.1  wavelengths,  the  mean  error  in  estimating  the  first-order  partial  is  -0.14 
dB  and  the  mean  emor  in  estimating  the  second-order  partial  is  -0.28  dB.  Note  that  at 
small  values  of  sensor  spacing,  the  standard  deviation  of  the  error  in  estimating  the 
second-order  partial  is  large  as  compared  to  the  mean  value  of  the  error.  The  opposite  is 
true  for  the  estimation  of  the  first-order  partial  derivative.  That  is,  the  standard  deviation 
of  the  error  is  small  as  compared  to  the  mean  value  of  the  error.  From  Figures  16  and  17, 
we  conclude  that  finite-differences  can  be  utilized  effectively  in  estimating  the  gradient 
of  the  pressure,  but  confidence  in  finite-difference  teehniques  deereases  significantly 
when  they  are  used  to  estimate  the  Hessian  of  the  pressure.  It  is  therefore  recommended 
that  an  acoustic  dyadic  sensor  be  realized  with  1 8  accelerometers  and  a  pressure  sensor  at 
the  origin.  This  will  eliminate  the  gradient  estimation  error  and  significantly  reduce  the 
Hessian  estimation  error,  since  the  Hessian  can  now  be  estimated  by  first-order  finite 
difference  approximations.  (Refer  to  equations  (14),  (15)  and  (16)). 
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Figure  16.  Statistical  behavior  of  the  error  in  estimating  the  mixed  second- 
order  partial  derivative  (with  respect  to  x  an  y)  at  the  origin. 
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Figure  17.  Statistics  of  the  error  in  estimating  the  first-order  (shown  at  the  left) 
and  pure  second-order  (shown  at  the  right)  partial  derivatives  of  the  pressure  at 
the  origin. 


5.0  Multichannel  Filtering  and  Directional  Acoustic  Sensors 

For  the  purposes  of  this  section,  it  is  convenient  to  deal  with  the  temporal  Fourier 
transform  of  the  Taylor  series  (8).  Generalizing  (8)  in  the  frequency  domain  we  obtain 

P{co,  T- )  =  ir,  (ro,  T )P{a),  P, )  +  IP\  {fo,  V )  ■  VP(fu,  i;, )  +  if,  {(o,  P )  ■  V VP(ru,  i;, )  •  tf,  (®,  F)  ,  (51) 

where  we  have  retained  only  the  first  three  terms  of  the  series  (dyadic  sensor)  and  the 
weights  and  W-,  arc  chosen  such  that  the  right  side  of  (5 1 )  is  consistent  with  the 

Taylor  series  approximation  of  the  pressure  field.  Specifically,  IF),  =  1 ,  W,  =  r  -r„  and 

I4P  =  “  '^1 )/ V2  are  the  only  weights  that  allow  (5 1 )  to  be  a  Taylor  series  extrapolation 

of  the  pressure  field  by  a  dyadic  sensor. 

It  is  important  to  note  that  (5 1 )  is  valid  for  any  analytic  acoustic  wave  field.  For 
the  special  case  of  an  arbitrary  plane  wave  we  have  that 

P(^o,P„)=  F(co)e'^  (52a) 

VP(co,  P„ )  =  JP(m,  F„  )/c  (52b) 

VVP(ro,  /), )  =  -P{co,  rj,  )kk  (52c) 


where  the  quantity  kk  is  a  dyadic  and  is  equivalent  to  the  rank  1  matrix  where 


A-  = 


k 


(53) 


Without  loss  of  generality,  we  can  let  r„  =  (0,0,0)  in  (52).  The  substitution  of  (52)  into 
(5 1 )  gives 


where 


P(fu,F)=F(ru,0^(A,F) 

(54) 

g(A',F)=l  +  7>-iy- 

(55) 

y  =  k  -  r  =  2k— Pi  •  r 

(56) 

The  magnitude  of  (55)  is 
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(57) 


4 


Equation  (54)  is  in  actuality  an  estimation  provided  by  a  dyadic  sensor  of  the  temporal 
Fourier  transform  of  the  pressure  of  a  plane  wave  at  the  field  point  specified  by  the  vector 
r  .  The  quality  of  the  estimation  would  depend  on  the  value  of  r/T  ,  where  /•  is  the 
distance  of  the  field  point  from  the  origin  and  T  is  the  wavelength.  Indeed,  from  (52)  and 
(54)  we  can  write 


fp  =20  log 


P{coJ-) 


=  201og|g(A'  -/-J 


(58) 


The  magnitude  of  the  g  function  is  seen  to  define  the  dB  error  resulting  from  estimating 
P{o}J-  )  by  means  of  replacing  the  Taylor  series  by  a  Taylor  polynomial  of  degree  2. 

Figure  18  contains  several  nonnalized  polar  plots  of  10''’^'"  =  pjp  plotted  as  a  function 
of  0^  for  (fij  =  (j)^  =  90'  and  0j  =  O'  where  the  angles  and  9^  appear  in  the  vector 


r  =  sin  cos  0^ x  +  sin  sin  6*5  j)  +  cos  z 


90  1  90  1 


90  1 


90  1 


Figure  1 8.  Polar  plots  of  the  normalized  Taylor  approximation  error  as  a 
function  of  Os  for  (|)s  =  (|)t  =  90^  and  Oj  =  0''  and  several  values  of  rlX, 


(59) 
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The  parameter  rj A  takes  on  values  of  2,  0.32,  0.2,  and  0.1.  The  plots  have  been 
normalized  so  that  the  maximum  value  is  unity.  Observe  that  the  polar  plots  in  Figure  18 
have  the  form  of  a  beam  pattern.  As  rj A  becomes  smaller,  the  beam  pattern  approaches 
an  omnidirectional  pattern.  The  r//i=2  pattern  has  a  beam  width  of  65.54"  Equation 
(58)  can  be  given  two  interpretations;  (1 )  Taylor  approximation  error  and  (2)  beam 
pattern.  Thus,  (54)  can  be  interpreted  as  the  generator  of  a  beam  pattern.  In  this 
intciprctation  the  value  of  /•  is  not  driven  by  the  need  for  an  acceptable  error  bound  on  the 
estimate  of  P[coJ-).  The  variable  r  simply  becomes  a  parameter  in  the  fomiing  of  a 

beam  and  the  function  g^(/r,r)  can  be  interpreted  as  a  discriminative  function.  The  vector 
f  in  (56)  becomes  a  steering  vector  and  we  are  now  interested  in  the  quality  of 
discrimination  of  the  beam  rather  than  the  quality  of  estimation  of  P{coJ-). 


From  Figure  1 8  we  sec  that  there  is  an  ambiguity  in  the  beam  patterns  (presence  of  two 
peaks).  Mathematically,  this  ambiguity  is  caused  by  the  fact  that  (57)  is  an  even  function 
of  Y .  It  would  be  beneficial  to  generate  beam  patterns  that  do  not  have  this  ambiguity. 
Beam  pattern  ambiguity  can  be  eliminated  or  significantly  reduced  by  picking  other  sets 
of  weights  in  (5 1 ).  However,  by  doing  so,  we  no  longer  have  a  valid  Taylor  series.  It  is 
more  correct  to  modify  (5 1 )  to 

5(m,  V)  =  r)P(co,  /";,)+  IF,  (ry,  f)  •  VP(co,  )  +  IF,  (gj,  F)-  VVP(co,  Ij, )  •  IF,  (co,  F)  (60) 


where  the  function  B{coJ-)  could  represent  the  output  of  a  frequency  domain  multi¬ 
channel  filter  with  P{coJ\^),  VP{(oJ\^)  and  WP{coJ\^)  as  the  inputs.  Refer  to  Figure  19 
for  a  pictorial  representation  of  this  concept. 


A'.r,) 

vvAri;,) 


>- 


Multi-channel  filter 


Figure  1 9.  Multi-channel  Filtering  and  Directional  Acoustic  Sensors 

Thus,  (60)  is  a  generalization  of  (51)  and  can  be  viewed  as  a  multi-channel  filtering 
approach  to  directional  sensors.  Consider  the  following  set  of  weights: 

W,{coJ-)=\ 

W,(£0,F)  =  a—-  (61) 

jeo  r 

W^{co,F)^^—- 
jeo  r 
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where  a  and  /?  are  two  parameters  that  may  be  a  function  of  frequency.  The  substitution 
of  (61)  into  (60)  results  in 

B{o},r)  =  p(cy,0^(«,  f)  (62) 

where 

g[fi,r)  =  \  +  a{h-f)+p{h-fY  (63) 

We  could  select  the  weights  in  (63)  to  steer  the  beam  of  the  directional  sensor  so  that  we 
have  maximum  sensitivity  in  the  look  direction  and  place  a  null  in  another  direction. 
Further,  we  could  select  the  weights  so  that  the  directional  sensor  produces  a  beam 
pattern  that  results  in  maximum  array  gain  against  isotropic  noise  [14].  Figure  20  shows 
the  two  weight-selection  criteria  for  the  vector  sensor.  The  maximum-array  gain  case  of 
Figure  20(a)  corresponds  to  a  =  3,yff  =  0 ,  whereas  the  optimum-null  case  of  Figure  20(b) 
corresponds  to  a  =  \,P  -0.  Figure  2 1  does  the  same  for  the  case  of  the  dyadic  sensor. 
The  maximum-array  gain  case  of  Figure  21(a)  corresponds  io  a  =  =  -5  ,  whereas 

the  optimum-null  case  of  Figure  21(b)  corresponds  to  a  =  2,P  =  \ .  Observe  that  in  each 
case,  the  ambiguity  was  eliminated  or  significantly  reduced. 


Figure  20.  Selection  of  weights  for  a  vector  sensor,  (a)  For 
maximum  array  gain,  (b)  For  optimum  null  placement. 


Array  Gain  «  9.5  dB  Array  Gain  -  7  dB 

(a)  (b) 

Figure  21.  Selection  of  weights  for  a  dya(iic  sensor,  (a) 
For  maximum  array  gain,  (b)  For  optimum  null 
placement. 
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6.0  Experimental  Results  Using  a  Partial  Dyadic  Sensor 

A  vector  sensor  placed  at  a  point  in  space  measures  both  the  pressure  and  the  gradient  of 
the  pressure  at  that  point.  A  dyadic  sensor  also  does  this.  In  addition,  it  also  measures  all 
of  the  spatial  second-order  partial  derivatives  of  the  pressure  at  the  point.  This  section 
will  discuss  a  partial  dyadic  sensor.  The  dyadic  sensor  is  partial  in  the  sense  that  it  only 
measures  the  spatial  pure  second-order  partial  derivatives  of  the  pressure  (i.e.,  the 
diagonal  elements  of  the  Hessian  matrix  (7)),  not  the  mixed  second-order  partials.  Figure 
5  shows  the  placement  of  the  components  of  the  partial  dyadic  sensor.  The  components 
consist  of  a  pressure  sensor  at  the  origin  and  six  CAVES-type  (flexural  disk) 
accelerometers  placed  at  the  extremities  of  the  three  orthogonal  lines  in  Figure  5. 

SITTEL  CORPORATION  used  the  services  of  EDO  Electro-Ceramic  Products,  Salt 
Lake  City,  Utah  to  build  the  partial  dyadic  sensor  illustrated  in  Figure  22.  Figure  23 
shows  the  support  structure  used  to  make  beam  pattern  measurements.  These 
measurements  were  performed  at  Seneca  Lake  in  May  2001. 

In  section  6. 1 ,  we  derive  the  theoretical  beam  patterns  for  the  elements  of  the  partial 
dyadic  sensor  and  provide  measured  results  from  Seneca  Lake  to  validate  the  theory.  In 
section  6.2,  we  show  how  multi-channel  signal  processing  of  a  vector  sensor  can  be  used 
to  realize  the  perfomiance  of  a  dyadic  sensor  when  the  pressure  field  is  assumed  to  be  an 
arbitrary  plane  wave.  We  also  provide  experimental  results  to  validate  this  claim. 


Figure  22.  Partial  dyadic  sensor  built  by  EDO  Electro-Ceramic  Products. 
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Figure  23.  Support  structure  used  for  partial  dyadic  sensor. 


6.1  Theoretical  and  Experimental  Beam  Patterns  for  a  Partial  Dyadic  Sensor 

The  goal  of  this  section  is  to  discuss  beam  patterns  produced  by  particular  accelerometer 
components  and  certain  combinations  of  accelerometer  components  of  the  EDO  sensor. 
Specifically,  we  will  focus  on  the  lines  labeled  p/dx^  and  d'^ p/dy'^  in  Figure  5.  The 
accelerometer  at  the  left  extremity  of  the  line  pjdx'  will  be  referred  to  as  the 
+  X  accelerometer,  whereas  the  accelerometer  at  the  right  extremity  will  be  called  the  -  jc 
accelerometer.  Similarly,  the  accelerometer  at  the  foremost  extremity  of  the  line 
d~ p I dy^  will  be  referred  to  as  the  +  y  accelerometer,  whereas  the  accelerometer  at  the 
rearmost  extremity  will  be  called  the  -  y  accelerometer.  Before  presenting  empirical 
beam  pattern  measurements,  it  is  necessary  to  first  discuss  the  theoretical  equivalents.  It 
is  well  known  in  the  theory  of  vector  sensors,  that  the  frequency-independent  beam 
pattern  of  a  single  accelerometer  follows  a  cos  0  {  +  x  or  -  x  accelerometer)  or  sin  0 
{  + y  or -y  aceelerometer)  pattern.  We  will  derive  the  theoretical  beam  pattern  for 

dajdx ,  which  can  be  approximated  by  the  first-order  finite  difference  of  the  -i-  x  and 
-  X  accelerometers.  Let  us  begin  with  the  relationship 


da^  _  I  d^p 
dx  p  dx~ 


(64) 
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This  relationship  is  valid  for  a  homogeneous  medium  resulting  in  a  constant  p .  For  a 
sinusoidal  plane  wave  of  frequency  (o ,  (64)  becomes 


da^tj')  k-  .  2  ,  2n 
— ^  =  — sm  cos  Bjp\t,r) 

dx  p 


The  mean-square  value  of  (65)  is 


where 


K  5x  j 


2  4  n 

COS  6j 


(l)jp-{t,r) 


(65) 


(66) 

(67) 


The  parameter  k  is  kept  constant  by  fixing  the  frequency  co  and  the  angle  (f)j.  We 
generate  a  theoretical  beam  pattern  for  da^/dx  by  plotting  (66)  as  a  function  of  the  angle 
0,  .  Similarly,  we  generate  a  theoretical  beam  pattern  for  dajdy  by  plotting 


2-4/1 

=  K  Sin  6j 


(68) 


as  a  function  of  the  angle  0j  . 

A  first-order  finite  difference  of  the  +  v  and  -  x  accelerometers  was  used  to  estimate  the 
left  side  of  (64).  Figures  24  (a)  and  (b)  give  the  empirical  and  theoretical  beam  patterns 
of  both  accelerometers  at  1200  Hz.  Both  empirical  patterns  are  close  to  the  theoretical 
“  cos^  ”  pattern.  Note  the  -25  dB  nulls  in  the  empirical  patterns.  Figure  24  (c)  consists  of 
the  empirical  and  theoretical  patterns  for  dajdx .  The  empirical  pattern  is  close  to  the 

theoretical  “  cos'  6  ”  pattern.  Note  the  -  2 1  dB  nulls  for  the  empirical  pattern  in  contrast 
to  the  infinite  nulls  for  the  theoretical  pattern.  Pulsing  a  source  at  a  fixed  frequency  of 
1200  Hz,  recording  and  processing  data,  and  then  rotating  the  vertical  shaft  in  Figure  23 
in  the  xy  plane  in  order  to  achieve  a  different  angle  6’^^  generated  the  patterns  for  da  jdx 

and  dajdy ,  respectively.  Figures  25  (a)  and  (b)  give  the  empirical  and  theoretical  beam 
patterns  of  the  -i-  y  and  -  y  accelerometers.  Both  empirical  patterns  are  close  to  the 
theoretical  “sin  0 ”  pattern.  Note  the  -25  dB  nulls  in  the  empirical  patterns.  Figure  25  (c) 
consists  of  the  empirical  and  theoretical  patterns  for  dajdy .  The  empirical  pattern  is 

close  to  the  theoretical  “siir  <9”  pattern.  Note  the  -  21  dB  nulls  for  the  empirical  pattern 
in  contrast  to  the  infinite  nulls  for  the  theoretical  pattern. 
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6.2  Achieving  Dyadic  Sensor  Performance  by  Multi-channel  Processing  of  a  Vector 
Sensor  in  a  Plane  Wave  Field 

In  this  section,  we  show  how  multi-channel  signal  processing  of  a  vector  sensor  can  be 
used  to  realize  the  performance  of  a  dyadic  sensor  when  the  pressure  field  is  assumed  to 
be  an  arbitrary  plane  wave.  We  also  provide  experimental  results  to  validate  this  claim. 

First,  consider  the  case  where  the  pressure  waveform  is  a  sinusoidal  (namowband)  plane 
wave  having  the  form 


/?(.v,  r,r,/)=  p{r,t)=  A +  ^ +  k  (69) 

The  gradient  and  the  Hessian  of  (69)  at  (v,, ,  y„ ,  2,, )  are 

V/;(,Y„ ,  v'„ ,z,^,t)  =  kA cosicot  +  (^  +  k  ■  i), )  (70) 

VV/;(,Y„ ,  v,| , 2„ , /)  =  -kk' A sin(r(;/  +  ^  +  k  -  p)  (71) 

where  /p  is  the  column  vector  of  wavenumbers.  Let  us  now  place  a  vector  sensor  at  the 
point  (x,| ,  v’l, ,  2|, ) .  This  sensor  measures 

Pi^'o  ^  .Vo ,  - 0  T )  =  H  sin (ry/  +  ^  +  A'  •  r„ )  (72) 

£(x„ ,  v„ , ^0 t)  =  — -kA cosfry/  +  i^  +  k  ■  )  (73) 

P 

where  the  column  vector  a_  contain  the  acceleration  components.  From  (73)  we  obtain 
the  3x3  matrix 


0  ^ 

£(T)  ,  3''o ,  -0  Tk'  (x,) ,  y,| ,  z„ ,  t)dt 


(74) 


Equations  (70)  and  (7 1 )  can  be  expressed  as 

Vpko  ^  To .  ^0  t)  =  -pa{x^ ,  y„ ,  ,  r) 


(75) 


VVjt?(Ao  ,  V() , j  ^2  ’  3^0  ’  ^0  ’  0  5  ^0  ?  (-^0  3  3^0  5  ^0  ?  0^^ 


The  above  pair  of  equations  indicates  that  for  time-haniionic  plane  wave  fields,  both  the 
pressure  gradient  and  the  pressure  Hessian  can  be  detennined  by  pressure  and 
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acceleration  measurements.  A  vector  sensor  located  at  (vf,, 
measurements. 

)  achieves  these 

Let  us  now  consider  an  arbitrary  (nan'owband  or  broadband)  plane  wave.  An  arbitrary 
plane  wave  moving  towards  the  origin  can  be  written  in  the  fomi 

/  t  f(  "-c 

p{x,}\Z,t)=f  t  + 

V  c  J 

(77) 

It  follows  Vp  and  S/Vp  have  the  values 

\  (  h-T‘\ 

Vp  =  n  g  t  + 

c  \  c  J 

(78) 

rvr-7  '  1  if  /?  •  F  ^ 

Wp~nn  —h  - 

C"  V  c  ) 

(79) 

where 

g  =  f 
h  =  g'  =  f" 

(80) 

The  Fourier  transforms  of  (77),  (78),  and  (79)  are 

F[pl=F’(7C<^)exp  JO) 

c 

(81) 

F[Vp]  =  « -yVt;  F[/?] 

c 

(82) 

F[vV;?]  =  -^'-!^6;-F[;?] 

(83) 

c 


The  Fourier  transform  of  the  acceleration  is  given  by 

F[£]  =  --fM  =  -— (84) 

P  pc 

Multiplying  (84)  by  the  transpose  of  its  conjugate  results  in 

=  (85) 

[pcy 
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Equation  (85)  leads  to  an  expression  for  nii  ,  namely. 


nn 


>1  =  {pc) 


2 


J  oy\¥\py  d  CO 


(86) 


The  quantities  V/;  and  VV/;  can  be  expressed  in  terms  of  the  pressure  and  acceleration  as 
follows: 


F[yp]  =  -pF\a] 


(87) 


F[vV/;]  = 


.  .  r  jF|£]F''bk^y 

p-co-F[p\~^  r  r - 

do) 


(88) 


Equations  (87)  and  (88)  indicate  that  for  an  arbitrary  plane  wave  field,  both  the  pressure 
gradient  and  the  pressure  Hessian  can  be  determined  by  pressure  and  acceleration 
measurements.  A  vector  sensor  located  at  (v,,,  )  achieves  these  measurements. 

Based  on  the  above  discussion,  we  can  use  equations  (60)  and  (61)  to  construct  a 
frequency  domain  multi-channel  filter  with  P{coJ\^),  VP{co,T\^)  and  VVP(<w,i;J  as  the 

inputs.  In  this  way,  we  can  obtain  the  performance  of  a  dyadic  sensor  with  vector 
sensor  m  easurem  ents. 

In  order  to  validate  this  claim,  tests  were  conducted  at  Seneca  Lake  during  the  period 
May  2001 .  A  vector  sensor  was  placed  at  the  origin  of  a  rectangular  coordinate  system 
and  a  400  Hz  acoustic  target  was  placed  atR,  =  =  90°relative  to  the 

vector  sensor  (refer  to  Figure  1 ).  Figure  26  shows  a  contour  plot  of  the  mean-squared 
value  of  (60)  as  a  function  of  frequency  and  azimuth  with  /?  =  0  (a  vector  sensor).  This  is 
the  Frequency-Azimuth  (FRAZ)  display  for  a  vector  sensor.  The  weights  of  the  multi¬ 
channel  filter  («  =  3,/?  =  0)  were  chosen  so  that  the  vector  sensor  beam  would  achieve 
maximum  array  gain  against  isotropic  noise.  Figure  27  shows  the  FRAZ  display  that 
results  when  the  weights  of  the  multi-channel  filter  {a  =  =  -5)  were  chosen  so  that 

the  equivalent  dyadic  sensor  beam  would  achieve  maximum  array  gain  against  isotropic 
noise.  A  comparison  of  Figures  26  and  27  shows  the  spatial  filtering  improvement  when 
using  the  extra  estimate  of  (88)  to  achieve  dyadic  perfonnance  from  a  vector  sensor. 

Using  the  same  test  set-up  as  above,  a  1300  Hz  interfering  source  was  placed  at 
R,  =  222  ft. ,9-1  =  225  ‘  =19'  relative  to  the  vector  sensor,  so  now  we  have  a  400 

Hz  target  and  a  1300  Hz  interferer  simultaneously  radiating  energy  toward  the  vector 
sensor.  Once  again,  a  comparison  of  Figures  28  and  29  shows  the  significant 
improvement  of  using  the  extra  estimate  of  (88)  to  aehieve  dyadic  perfonnance  from  a 
vector  sensor. 
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Figure  26.  FRAZ  display  using  vector  sensor  processing  (single  400  Hz  target 
located  at  r  =  130  ft,  0t  =  0”,  =  90”). 
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Figure  27.  FRAZ  display  using  dyadic  sensor  processing  (single  400  Hz  target 
located  at  r  =  130  ft,  Oj  =  0°,  ^  =  90”). 
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Figure  28.  FRAZ  display  using  vector  sensor  processing  (400  Hz  target  located  at 
r  =  130  ft,  O'l'  =  0",  (|)t  =  90”and  1300  Hz  interfercr  located  at  ri  =  222  ft,  0|  =  225^ 
^1  =  19^^). 
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Figure  29.  FRAZ  display  using  dyadic  sensor  processing  (400  Hz  target  located  at  r  =  130  ft, 
Oj  =  0^  (j)!'  =  90",  and  1300  Hz  interferer  located  at  V]  =  222  ft,  0i  =  225",  ({>1  ^  )• 
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7.0  Effects  of  Quantization  and  Wide  Band  Noise  on  Finite  Differences  Between 
Two  Accelerometers 


This  section  analyzes  the  effeets  of  quantization  and  wide  band  noise  on  the  finite 
differences  between  sinusoidal  signals  received  from  two  accelerometers  spaced  three  to 
six  inches  apart.  The  noise  bandwidth  is  2000  Hz  and  is  comprised  of  the  following  three 
components: 


(1 )  Ambient  Ocean  noise 

(2)  Preamplifier  noise 

(3)  Analog  to  Digital  Converter  (ADC)  quantizing  noise. 


Ambient  noise  ranges  from  the  maximum  acceptable  noise  as  shown  in  the  NUWC 
Vector  Sensor  Specifications  dated  03  Feb  2000,  to  a  value  40dB  higher.  An 
accelerometer  sensitivity  of -18  dBv/g  is  used.  Preamplifier  noise  is  -166dBv/VHz, 
which  is  equivalent  to  the  noise  floor  of  a  SITTEL  CORPORATION  signal  conditioner. 
ADC  quantizing  noise  is  lO.SdB  below  the  least  significant  bit  (LSB)  level  or  quantizing 


Accelerometer  1  Gain:  40dB  16  to  24  Bits 

Sens: -18dBv/g  BW:  50-2050Hz  1  Volt  Full  Scale 


Accelerometer  2  Gain:  40dB  16  to  24  Bits 

Sens:-18dBv/g  BW:  50-2050Hz  1  Volt  Full  Scale 


Finite 

Difference 


Figure  30.  Modeling  Setup 


step  size  of  the  ADC.  For  an  N-bit,  bipolar  ADC  with  full-scale  amplitude  of  ±  Fs,  the 
LSB  level  is  Fs/2*^''’.  The  analysis  uses  both  a  16-bit  ADC  and  24-bit  ADC  to  observe 
the  effects  of  quantizing  step  size  on  the  finite  difference  errors.  Figure  30  shows  the 
modeling  setup  used  for  the  analysis. 

The  analysis  is  comprised  of  the  following  four  parts: 
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Part  1.  The  signal  frequency  is  set  at  50Hz,  which  is  considered  to  be  worst  case 
since  the  finite  differences  are  the  smallest  at  the  lowest  signal  frequencies  within  the 
band  with  fixed  accelerometer  spacing.  Noise  is  considered  to  be  uncorrelated  between 
accelerometers.  Accelerometer  spacing  is  3  inches  and  6  inches.  The  RMS  value  of  the 
composite  noise  ranges  from  1  LSB  to  100  LSB’s  for  the  16-bit  ADC  and  from  256 
LSB’s  to  25,600  LSB’s  for  the  24-bit  ADC  with  the  sinusoidal  signal  levels,  in  each  case, 
ranging  from  OdB  to  +50dB  above  the  composite  noise  level.  The  signals  are  assumed  to 
be  plane  wave  and  arrive  along  the  axis  formed  by  the  two  accelerometers. 

Part  2.  The  signal  frequency  is  set  at  50Hz.  Noise  is  uncorrelated  between 
accelerometers.  Accelerometer  spacing  is  3  inches  and  6  inches.  The  RMS  value  of  the 
composite  noise  ranges  from  1  LSB  to  100  LSB’s  for  the  16-bit  ADC  and  from  256 
LSB’s  to  25,600  LSB’s  for  the  24-bit  ADC  with  the  sinusoidal  signal  levels,  in  each  case, 
ranging  from  OdB  to  +50dB  above  the  composite  noise  level.  The  signals  are  assumed  to 
be  plane  wave  and  arrive  along  the  axis  formed  by  the  two  accelerometers.  The  signal  is 
sensed  with  2  accelerometer  sets  and  averaged  to  determine  improvement  in  errors  due  to 
averaging. 

Part  3.  The  signal  frequency  is  set  at  50Hz.  The  noise  correlation  coefficient 
between  accelerometers  ranges  from  0%  to  95%  to  determine  the  improvements  in  finite 
difference  errors  due  to  correlated  noise.  The  noise  and  signal  attributes  are  as  described 
in  part  1  and  part  2.  A  16-bit  ADC  converter  is  used  with  3-inch  spacing  between 
accelerometers.  With  close  accelerometer  spacing,  ambient  noise  is  likely  to  be 
correlated  as  described  in  Appendix  A. 

Part  4.  The  signal  frequency  varies  from  50Hz  to  2000Hz  with  a  constant 
accelerometer  spacing  of  3  inches.  Improvement  is  expected  at  higher  frequencies  due  to 
greater  finite  differences.  The  ADC  is  set  at  16  Bits  with  the  input  noise  level  set  at  lOdB 
above  the  LSB.  The  signal-to-noise  ratio  varies  from  OdB  to  50dB. 

Conclusions.  It  is  feasible  to  measure  finite  differences  of  50Hz  sinusoidal  signals,  in 
the  presence  of  wide  band  noise,  emanating  from  accelerometers  with  3-inch  spacing 
between  accelerometers.  When  input  signal  plus  noise  levels  are  large  enough  to  dither 
the  ADC,  the  analysis  shows  that  the  errors  depend  upon  the  signal  to  noise  ratio.  An 
acceptable  dithering  level  is  lOdB  or  greater  relative  to  the  LSB  of  the  ADC.  At  this 
level,  there  is  no  difference  in  error  perfonnance  of  a  16-bit  ADC  and  a  24-bit  ADC. 
Normalized  eirors  are  the  difference  between  the  perturbed  finite  difference  and  actual 
finite  difference  divided  by  the  actual  finite  difference. 

For  a  50Llz  signal  from  2  accelerometers  spaced  3  inches  apart,  with  wide  band  noise 
uncorrelated  between  accelerometers,  a  SNR  of  30dB  is  required  to  obtain  a  normalized 
average  RMS  error  of  10%.  The  normalized  average  RMS  error  is  the  square  root  of  the 
normalized  mean  square  error  averaged  over  a  cycle  of  signal  frequency.  For  a 
normalized  maximum  RMS  error  of  10%,  the  SNR  must  be  in  excess  of  50dB.  The 
normalized  maximum  RMS  error  is  the  square  root  of  the  maximum  normalized  MSB 
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within  a  cycle  of  signal  frequency.  In  this  analysis,  the  maximum  MSE  occurs  at  the  peak 
of  the  sinusoidal  signal  where  the  minimum  finite  difference  occurs.  The  minimum  finite 
difference  is  defined  as  A*(l-cos  (A0))  where  A  is  the  peak  signal  amplitude  and  A0  is 
equal  to  2*n:*d/X,,  with  d  equal  to  the  distance  between  accelerometers  and  X  equal  to  the 
wavelength.  A  SNR  of  36dB  is  required  for  a  normalized  maximum  RMS  eiTor  equal  to 
the  minimum  finite  difference. 

The  analysis  has  shown  that  increasing  the  spacing  between  the  accelerometers  results  in 
lower  errors  for  the  same  signal  frequency.  At  50Hz,  a  6dB  improvement  is  realized  by 
6-inch  spacing  over  3-inch  spacing. 

The  analysis  has  shown  that  averaging  the  finite  differences  of  2  sets  of  accelerometers 
looking  at  the  same  signal  results  in  3dB  error  improvement  as  long  as  the  noise  is 
uncorrelated  between  sets  of  accelerometers. 

The  analysis  has  shown  that  lower  eiTors  occur  when  the  noise  between  accelerometers  is 
correlated.  The  improvement  is  7dB  from  0%  correlation  to  90%  correlation  at  a  given 
SNR.  The  actual  noise  cori'elation  is  not  known  but  is  expected  to  be  high  due  to  the 
close  spacing  of  the  accelerometers. 

Finally,  the  analysis  has  shown  that  50Hz  has  the  worst  case  errors.  The  errors  for  3-inch 
spacing  improve  at  higher  frequencies.  For  frequencies  up  to  50OHz,  the  error 
improvement  is  6dB  per  octave  of  frequency.  For  frequencies  above  SOOHz,  the 
improvement  is  substantial,  although  less  than  6dB  per  octave. 

Analysis.  A  Monte  Carlo  analysis  was  used  to  compute  the  difference  error  due  to  noise 
and  quantizing  effects.  A  signal  cycle  was  divided  into  400  phases  with  each  phase 
equivalent  to  the  phase  shift  between  accelerometers.  The  finite  difference  is  the 
difference  in  instantaneous  amplitudes  between  2  adjacent  phases  of  the  signal.  The 
mean  square  error  (MSE),  between  the  finite  difference  corrupted  by  noise  and  ADC 
effects,  and  the  actual  finite  difference,  was  computed  at  each  of  the  400  phases  using 
2048  sample  functions  of  noise  at  each  ADC  input.  The  average  MSE  was  then 
computed  by  averaging  the  MSE  over  a  cycle  of  signal  input  and  normalizing  to  the 
actual  finite  difference.  Likewise  the  maximum  MSE  was  computed  by  determining  the 
maximum  MSE  over  a  cycle  and  normalizing  to  the  actual  finite  difference.  It  is 
expected  that  the  maximum  error  occurred  at  the  minimum  finite  difference.  As  shown 
in  Appendix  A,  the  minimum  finite  difference  is  equal  to  A*(l-cos  (A0)),  where  A  is  the 
peak  signal  amplitude  and  A0  is  equal  to  2*n*d/X,  with  d  equal  to  the  distance  between 
accelerometers  and  X  equal  to  the  wavelength.  The  maximum  finite  difference  is  equal  to 
A*sin  (A0). 

With  3  inch  spacing,  at  50Hz,  A0  is  equal  to  0.9  degrees,  and  the  minimum  finite 
difference  is  equal  to  A*  123*  10'°^.  The  ratio  of  maximum  to  minimum  finite  difference 
is  127  or  42.1  dB.  With  6  inch  spacing,  at  50Hz,  A0  is  equal  to  1.8  degrees,  and  the 
minimum  finite  difference  is  equal  to  A*493*10'^^.  The  ratio  of  maximum  to  minimum 
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finite  difference  is  63  or  36.1  dB.  The  finite  differences  are  directly  proportional  to 
frequency  for  a  constant  spacing  between  accelerometers. 


Analysis  Part  1 

Figures  3  1  and  32  show  the  normalized  average  MSE  and  maximum  MSE,  respectively, 
versus  input  noise  level  and  signal  to  noise  ratio  (SNR)  for  a  16  bit  ADC  and  3-inch 
spacing.  The  minimum  noise  level  is  about  lOdB  above  the  ADC  LSB.  A  30dB  SNR  is 
required  for  a  normalized  MSE  of  1%,  which  corresponds  to  a  normalized  RMS  eiror  of 
10%.  However,  for  a  normalized  maximum  MSE  of  ]%,  a  SNR  in  excess  of  50dB  is 
required.  A  36dB  SNR  is  required  for  a  normalized  maximum  MSE  less  than  the 
minimum  finite  difference. 

These  plots  show  that  the  errors  are  constant  with  input  noise  levels  and  a  function  of 
only  the  SNR  for  input  levels  that  are  lOdB  or  more  above  the  LSB  of  the  ADC.  The 
slight  decrease  of  errors  at  levels  below  lOdB  relative  to  the  LSB,  for  SNR’s  OdB  to 
30dB  arc  unexplained  and  may  be  anomalies  caused  by  small  input  levels  relative  to  the 
ADC  LSB.  Since  it  is  good  practice  to  maintain  minimum  RMS  noise  into  an  ADC  at  a 
level  lOdB  or  more  above  the  LSB,  the  data  below  lOdB  is  ignored. 

Figures  33  and  34  show  the  normalized  average  MSE  and  maximum  MSE,  respectively, 
versus  input  noise  level  and  SNR  for  a  24  bit  ADC  and  3-inch  spacing.  The  minimum 
noise  level  is  about  48dB  above  the  ADC  LSB.  It  can  been  seen  that  the  errors  are 
equivalent  to  those  of  the  16  bit  ADC  with  input  noise  13dB  or  higher  relative  to  the 
LSB. 

It  can  be  concluded  that  as  long  as  the  input  RMS  noise  level  is  kept  above  the  LSB  of 
the  ADC,  the  finite  difference  errors  of  a  16  bit  ADC  is  comparable  to  those  of  a  24  bit 
ADC.  A  30dB  SNR  is  required  to  keep  the  average  MSE  below  10%.  A  36dB  SNR  will 
keep  the  maximum  MSE  below  the  minimum  finite  difference. 

Figures  35  and  36  show  the  average  MSE  and  maximum  MSE,  respectively,  for  a  spacing 
of  6  inches,  with  50Flz  input  into  a  16-bit  ADC.  As  shown  in  the  figures,  the  errors  are 
6dB  lower  than  those  with  3-inch  accelerometer  spacing  are.  This  suggests  that  higher 
frequencies  will  have  lower  errors  at  3-inch  spacing.  Errors  as  a  function  of  frequency 
are  addressed  in  part  4  of  this  analysis. 

Analysis  Part  2 

Figures  37  and  38  show  the  average  MSE  and  maximum  MSE,  respeetively,  for  a  spacing 
of  3  inches,  with  50Hz  input  into  a  16-bit  ADC  with  the  finite  differences  obtained  from 
an  average  of  2  accelerometer  sets.  As  shown  in  the  figures,  the  errors  are  3dB  lower 
than  those  obtained  with  a  single  set  of  accelerometers.  This  analysis  shows  that  as  long 
as  the  noise  is  uncorrelated  between  sets  of  accelerometers,  an  improvement  in  finite 
difference  eiror  can  be  realized  by  averaging. 
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Analysis  Part  3 


Adding  portions  of  uncorrelated  noise  from  one  accelerometer  into  the  other 
accelerometer  simulated  correlated  noise.  Defining  N1  as  the  noise  from  accelerometer 
#1  and  N2  as  the  noise  from  accelerometer  #2,  the  correlated  noise  inputs  are  defined  as 
NT  =  N1  +  a*N2,  and  N2’  =  a*Nl  +  N2,  where  a  is  a  coefficient,  less  than  unity 
defined  by  the  correlation  coefficient,  p,  as  follows:  a  =  [l-(l-p  )  ]/p.  Figures  10  and 
11  show  the  improvements  in  average  MSB  and  maximum  MSB,  respectively,  due  to 
noise  correlation.  For  noise  correlated  at  90%,  a  7dB  improvement  is  realized. 

It  can  be  concluded  that  noise  correlation  will  result  in  an  improvement  in  finite 
differences. 

Analysis  Part  4 

The  signal  frequency  was  varied  from  50Hz  to  2000Hz  with  a  constant  accelerometer 
spacing  of  3  inches.  The  ADC  was  set  at  16  Bits  with  the  input  noise  level  set  at  lOdB 
above  the  LSB.  The  signal-to-noise  ratio  was  varied  from  OdB  to  50dB.  Figures  41  and 
42  show  the  improvements  in  average  MSB  and  maximum  MSB,  respectively,  as  a 
function  of  frequency.  At  lower  frequencies  (f  <  500),  the  improvement  is  6dB/octave  of 
frequency.  At  higher  frequencies,  the  improvement  is  less  than  6dB/octave. 

It  can  be  concluded  that  50Hz  signals  in  the  presence  of  uncorrelated  wide  band  noise  is 
the  worst  case  for  finite  difference  errors. 

Recommendations.  This  analysis  suggests  that  a  dyadic  sensor,  covering  the  frequency 
band  of  50  Hz  to  2,050  Hz,  can  be  built  in  a  small  package  like  the  form  factor  of  the 
existing  DT  276  hydrophone  (a  3-inch  diameter  by  6-inch  long  right  circular  cylinder). 
In  order  to  fit  the  necessary  accelerometers  in  this  form  factor,  SITTBB  recommends  the 
use  of  singlecrystal  or  MBMS  sensor  technology.  SITTBB  is  currently  developing  a 
complete  dyadic  sensor  and  signal  processing  design  using  this  advanced  sensor 
technology. 
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Figure  37.  Average  MSB  with  16  Bit  ADC, 
3-Inch  Accelerometer  Spacing,  and  Average 
of  2  Accelerometer  Sets 
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Figure  39,  Average  MSB  with  16  Bit  ADC 
as  a  Function  of  Noise  Correlation 
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Figure  38.  Maximum  MSB  with  16  Bit  ADC, 
3-Inch  Accelerometer  Spacing,  and  Average 
of  2  Accelerometer  Sets 
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Figure  40.  Maximum  MSB  with  16  Bit  ADC 
as  a  Function  of  Noise  Correlation 
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Figure  41 .  Average  MSB  with  16  Bit  ADC 
and  3  Inch  Accelerometer  Spacing,  as  a 
Function  of  Frequency 
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Figure  42.  Maximum  MSB  with  24  Bit  ADC 
and  3  Inch  Accelerometer  Spacing,  as  a 
Function  of  Frequency 
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8.0  Summary  and  Recommendations  For  Future  Research 
8.1  Summary 

In  Section  2.0,  we  formulated  a  rigorous  mathematical  framework  for  the  general 
theory  of  directional  acoustic  sensors.  Specifically,  the  number  of  terms  in  the  multi¬ 
dimensional  Taylor  series  expansion  of  the  acoustic  pressure  field  about  some 

measurement  point  Fg  defines  the  order  of  the  directional  sensor.  A  scalar  acoustic 
pressure  sensor  (e.g.,  a  hydrophone)  was  defined  as  a  directional  acoustic  sensor  of  order 
zero.  This  sensor  only  measures  the  scalar  acoustic  pressure  p{t,ro)  at  the  point  Fp ;  its 
Taylor  series  about  /;,  assumes  that  the  acoustic  pressure  field  p{t,r)  about  that  point  is 
independent  of  the  field  point  F  =  (x,y,z).  An  acoustic  vector  sensor  was  defined  as  a 
directional  acoustic  sensor  of  order  one.  This  sensor  measures  both  /?(t,Fo)  and  the 
pressure  gradient/vector  V/;(t,Fo)  at  the  point  F^, ;  its  Taylor  series  about  Fq  assumes  that 
the  acoustic  pressure  field  p{t,r)  about  that  point  is  a  linear  function  of  the  field  point  F . 
Similarly,  an  acoustic  dyadic  sensor  was  defined  as  a  directional  acoustic  sensor  of 
order  two.  This  sensor  measures  pit,?,^),  Vp(t,Fo)  and  the  dyadic  Wp{t,rQ)  at  the  point 
z;, ;  its  Taylor  series  about  assumes  that  the  acoustic  pressure  field  p{t,r)  about  that 

point  is  a  quadratic  function  of  the  field  point  F . 

Equation  (8)  is  the  pressure  estimate  generated  by  the  dyadic  sensor  and  is  capable  of 
extrapolating  the  acoustic  pressure  field  beyond  the  measurement  point  so  that  it 

actually  knows  this  field  at  every  point  inside  a  sphere  of  radius  ./?  =  |F  -  zi| ,  where  R  is 

defined  by  the  type  of  directional  sensor  and  the  error  associated  with  the  extrapolation. 
The  Taylor  polynomial  (8)  is  the  formula  used  to  do  the  wave  field  extrapolation  with 
some  specified  error  The  normalized  mean-squared  value  of  s(t,r)  vs.  H/ A  is 

plotted  in  Figure  2.  For  a  10%  pressure  estimation  error,  the  dyadic  sensor  measures  the 
pressure  field  everywhere  inside  a  sphere  of  radius  R  =  A/ 4 .  In  this  sense,  the  dyadic 
sensor  is  equivalent  to  a  volumetric  spherical  array  of  pressure  sensors.  In  practical 
terms,  this  means  that  for  an  acceptable  pressure  estimation  error,  a  dense  volumetric 
array  of  pressure  sensors  can  be  replaced  by  a  sparse  array  of  dyadic  sensors. 

In  Section  3.0,  we  answered  the  question:  If  three  accelerometers  and  a  pressure 
sensor  (i.e.,  a  vector  sensor)  can  achieve  a  maximum  array  gain  against  isotropic  noise 
of  6.0  dB,  a  question  that  naturally  follows  is  how  many  more  individual  sensors  (e.g., 
accelerometers  and/or  pressure  sensors)  must  fit  into  a  given  sensor  housing  or  form 
factor  to  improve  upon  the  aforementioned  directionality  of  the  acoustic  vector  sensor? 
Recall  that  the  dyadic  sensor  measures  p{t,f),  V/>(i,Fj  and  the  dyadic  VV;7(^,Fo)  at  the 
point  .  According  to  (6),  (7),  and  (8),  we  must  measure  13  terms  involving  the  pressure 

and  the  first  and  second  partial  derivatives  of  the  pressure.  Due  to  the  symmetry  of  the 
mixed  partials  in  (7),  the  number  of  distinct  terms  reduces  to  10.  We  demonstrated  that 
these  10  terms  can  be  estimated  by  1 9  judiciously  spaced  pressure  sensors,  or  18 
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judiciously  spaced  accelerometers  and  1  pressure  sensor  placed  at  the  origin.  Under  these 
conditions,  the  dyadic  sensor  could  achieve  a  maximum  array  gain  against  isotropic  noise 
of  9.5  dB. 

In  Section  4.0,  we  examined  the  effects  of  finite-difference  approximations  on  the 
estimation  of  the  pressure  gradient  and  Hessian  (realization  by  pressure  sensors)  and  also 
examined  the  effects  of  f  nite-difference  approximations  on  the  estimation  of  Va 
(realization  by  accelerometers).  Specifically,  this  section  analyzed  some  of  the  errors  that 
occur  when  finite-differences  are  used  to  estimate  the  first-order  and  second-order  partial 
derivatives  of  the  aeoustic  pressure.  The  errors  considered  were  those  caused  by  a  (1) 
finite,  nonzero  spacing  A.v  between  sensors  and  (2)  an  amplitude  mismatch  among  the 
composite  frequency  responses  of  the  various  sensors.  Expressions  for  these  errors  were 
derived  and  plots  were  generated  that  illustrated  the  dependency  of  error-behavior  on 
relative  sensor  spacing  and  dB-mismatch.  From  Figures  16  and  17,  we  conclude  that 
finite-differences  can  be  utilized  effectively  in  estimating  the  gradient  of  the  pressure,  but 
confidence  in  finite-difference  techniques  decreases  significantly  when  they  are  used  to 
estimate  the  Hessian  of  the  pressure.  It  is  therefore  recommended  that  an  acoustic  dyadic 
sensor  be  realized  with  18  accelerometers  and  a  pressure  sensor  at  the  origin.  This  will 
eliminate  the  gradient  estimation  error  and  significantly  reduce  the  Hessian  estimation 
error,  since  the  Hessian  can  now  be  estimated  by  first-order  finite  difference 
approximations.  (Refer  to  equations  (14),  (15)  and  (16)). 

Section  5.0  discussed  how  the  dyadic  sensor  can  be  viewed  as  a  multi-channel  filter. 
Refer  to  Figure  19.  Specifically,  equation  (60)  is  the  temporal  Fourier  transform  of  the 
output  of  the  filter  and  the  inputs  to  the  filter  are  the  temporal  Fourier  transforms  of  the 

pressure,  pressure  gradient  and  Hessian  of  the  pressure.  The  filter  weights  Wf,,  W^  and 

can  be  functions  of  frequency  and  can  be  chosen  so  that  the  filter  output  is  either  an 
estimation  of  the  pressure  (a  Taylor  polynomial),  a  beam  that  achieves  maximum  array 
gain  against  isotropic  noise,  or  a  beam  that  can  be  shaped  in  some  desired  fashion,  like 
placing  a  null  in  the  beam  pattern  to  spatially  filter  out  an  interfering  source.  With  proper 
selection  of  the  weights  to  achieve  maximum  array  gain  against  isotropic  noise,  a  dyadic 
sensor  can  produce  a  65-degree  beam  (3-dB  beam  width).  Under  the  same  conditions,  a 
vector  sensor  can  produce  a  1 05-degree  beam. 

The  experimental  results  in  Section  6. 1  show  that  the  elements  of  a  partial  dyadic 
sensor  can  produce  acceptable  beam  patterns  with  nulls  between  -20  to  -25  dB.  The 
experimental  results  in  Section  6.2  show  that  a  dyadic  sensor  is  a  very  effective  multi¬ 
channel  spatial  filter. 

Section  7.0  showed  that  first-order  finite-difference  approximations  can  be 
successfully  performed  in  the  presence  of  quantization  and  wide  band  noise  effects.  A 
worst-case  analysis  shows  that  using  two  accelerometers  spaced  3  inches  apart,  one  can 
achieve  a  normalized  RMS  error  of  10%  when  attempting  a  first-order  finite  difference  of 
a  50  Hz  plane  wave  in  a  2000  Hz  band.  This  requires  a  signal-to-noise  ratio  (SNR)  of  at 
least  30  dB  and  assumes  that  the  noise  at  each  accelerometer  is  uneorrelated.  This  can  be 
done  using  a  16-bit  analog-to-digital  converter  (ADC).  Further,  for  frequencies  up  to 
2000  Hz,  the  required  SNR  drops  off  at  about  6  dB  per  octave.  This  implies  that  the 
same  3-inch  spaced  accelerometers  only  need  a  SNR  of  6  dB  to  achieve  the  same  10  % 
error  for  a  first-order  finite  difference  of  an  800Hz  plane  wave. 
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8.2  Recommendations  for  Future  Research 


Based  on  the  results  in  this  report,  SITTEL  CORPORATION  recommends  the 
design,  build  and  test  of  two  complete  dyadic  sensors;  each  sensor  using  18  underwater 
accelerometers  and  a  single  pressure  sensor  housed  in  the  existing  DT  276  form  factor. 

In  order  to  minimize  the  finite  difference  approximation  errors  and  fit  within  the  DT  276 
form  factor,  SITTEL  recommends  the  use  of  single  crystal  or  MEMS  technology  to 
implement  the  accelerometers.  SITTEL  intends  to  perform  extensive  element  pattern 
measurements  associated  with  the  complete  Hessian  matrix  and  compare  these 
measurements  to  theoretical.  For  a  general  linear  acoustic  wave,  we  will  perform 
extensive  signal  processing  algorithms  to  show  the  improvement  over  the  FYOl  partial 
dyadic  sensor  and  vector  sensor.  We  will  also  demonstrate  how  two  dyadic  sensors  can 
improve  an  intensity-based  ranging  algorithm  using  vector  sensors.  Further,  we  will 
show  that  for  a  single  plane  wave  source,  the  outputs  of  a  dyadic  sensor  can  be  processed 
in  a  triadic  signal-processing  algorithm  to  achieve  triadic  sensor  perfomiance. 
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APPENDIX  A 


1.  Finite  Differences. 

Define:  Si  =  A*sin  (cot) 

S:  =  A*sin  (03t  -  kd),  where:  k  =  2*kIX  and  d  =  spacing  between 

accelerometers 

Let:  cot  =  0  and  kd  =  A0 

Then:  S|  -  S2  =  A*(sin0  -  sin  (0  -  A0))  which  is  maximum  at  0  =  0 

And  0  at  0  =  (71  -  0)/2 

The  minimum  step  size  used  in  this  analysis  occurs  when  0  =  nil. 

In  this  case  S|  -  S2  =  A*(l -cos  (A0)) 

2.  Effects  of  correlation. 

Define:  X|  =  S|+N| 

X2  =  S2  +  N2 

With:  a.s"  =  signal  power,  equal  at  both  accelerometers 

CT„"  =  noise  power,  equal  at  both  accelerometers 
Ps  =  signal  correlation  coefficient 
p„  =  noise  correlation  coefficient 


X|  -X2  =  S|  -S2  +  N,  -N2 

Output  power  =  <(X  1  -  X2)“>  =  <(S  1  -  S2)'>  +  <(  N 1  -  N2)">, 

=  2as‘(l  -ps)  +  2an^(l  -  p,,) 

Where  <>  denotes  expected  value.  Signal  cross  noise  terms  are  equal  to  0  since  signal  is 
uncorrrelated  to  noise. 

If  signal  is  sinusoidal  ps=  cos  (kd)  or  if  broadband  over  range  F|  to  F2, 

Ps  =  cos  (cooT)*sin  (7iBx)/(7rBT),  with  coo  =  ti*(F2+F|),  x  =  d/c  and 
B  =  F2-F| 

If  noise  is  isotropic  p„  =  sin  (kd)/kd, 

For  broadband  noise: 

Ay/ 2 

<pn>  =  (kd2-kd|)'’*  \d(kcl)* Sin{kd)l(kd) ,  for  d  =  3”,  F|  =  50,  F2  =  2050, 

Ay/1 

<p„>  =  0.976  and  2an^(l  -  <p„>)/2an^  =  -16dB  or  noise  is  reduced  by  16dB  due  to 
correlation 
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